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ABSTRACT 

The distribution of early-type galaxy velocity dispersions, </>(cr), is measured using a sample drawn 
from the SDSS database. Its shape differs significantly from that which one obtains by simply using the 
mean correlation between luminosity L and velocity dispersion a to transform the luminosity function 
into a velocity function: ignoring the scatter around the mean a — L relation is a bad approximation. An 
estimate of the contribution from late-type galaxies is also made, which suggests that 4>(cr) is dominated 
by early-type galaxies at velocities larger than ~ 200 km/s. 

Subject headings: galaxies: elliptical — galaxies: evolution — galaxies: fundamental parameters — 
galaxies: photometry — galaxies: stellar content 



1. INTRODUCTION 

The distribution of internal velocities in galaxies, the 
"velocity dispersion function" , figures prominently in sev- 
eral cosmological calculations. Gravitational lensing cross 
sections (Turner, Ostriker and Gott 1984) and the masses 
of central black holes (Ferrarese & Merritt 2000; Gebhart 
et al. 2000) depend more closely on the internal velocities 
than on total mass or total light. The lensing efficiency of 
a galaxy is expected to scale as the fourth power of its ve- 
locity dispersion, and recent work suggests that the mass 
of the central black hole scales as the fourth or fifth power 
of the velocity dispersion of the host galaxy (Tremaine et 
al. 2002). Therefore, a reliable estimate of the velocity 
function is extremely useful. 

Perhaps more importantly, the velocity function pro- 
vides a crucial link between models for galaxy formation 
and the observed universe. Although the shape of the mass 
function of dark matter halos is routinely predicted by the- 
oretical models (Press & Schechter 1974; Sheth & Tormen 
1999; Jenkins et al. 2001), a comparison with the mass 
function of galaxies is not straightforward. This is because 
sufficiently massive halos host more than one galaxy, so 
there is no simple correspondence between a halo's mass 
and the masses of the galaxies it hosts. Although galaxy 
formation models (White & Rees 1978; White & Frenk 
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1991) predict the numbers and masses of galaxies which 
form in halos as a function of parent halo mass (e.g., Kauff- 
mann et al. 1999; Somerville & Primack 1999; Springel et 
al. 2001; Benson et al. 2002), so they can be used to 
predict the mass function of galaxies, the total mass of 
a galaxy is notoriously difficult to measure (but see, e.g., 
Kauffmann et al. 2003 for a technique which estimates 
the contribution to the mass which comes from stars). 
The same theoretical models also predict the distribution 
of luminosities, the galaxy luminosity function, which is 
much easier to measure (see, e.g., the recent determina- 
tions by Cross et al. 2001; Blanton et al. 2002; Madgwick 
et al. 2002). However, this comparison between theory 
and observation depends, of course, on waveband. The 
distribution of internal velocity dispersions (for early-type 
galaxies) , or of circular velocities (for spiral galaxies) does 
not depend on waveband. Since the theoretical models do 
predict </>(c), the shape of the velocity function of galax- 
ies, it is a more direct way of testing models of galaxy 
formation than is the luminosity function <j)(L). 

Shortly after the first measurements of velocity disper- 
sions in other galaxies were made (e.g., Minkowski 1954), 
it was recognized that luminosity and velocity disper- 
sion of early- type galaxies are correlated (e.g., Poveda 
1961; Minkowski 1962; Fish 1964; Faber & Jackson 1976): 
a cx L 1 ^ , with the exact value of ip depending on wave- 
length. And, following Tully & Fisher (1977), a simi- 
lar correlation between the luminosity and circular veloc- 
ity v c of spiral galaxies has also been extensively stud- 
ied (e.g., Giovanelli, Hayncs et al. 1997; Verheijen 2000). 
Some (e.g., Shimasaku 1993; Gonzalez et al. 2000) have 
used these correlations to convert a measured distribu- 
tion of luminosities <p(L) into an estimate of 4>(v c ), sim- 
ply by using the v c — L relation to transform variables: 
4>( v c) = <fi{L)\dL/dv c \. This procedure ignores the fact 
that there is scatter in the v c — L or a — L relations, and so 
it almost certainly underestimates the number of objects 
which have large velocity dispersions. Since objects with 
large velocity dispersions figure prominently in a number 
of different arguments regarding gravitational lensing and 
galaxy formation, simply changing variables is an unreli- 
able way to proceed (e.g., Kochanek 1993). 
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The majority of lenses are known to be early-type galax- 
ies. Recently, Bernardi et al. (2003a, b,c,d) compiled 
a sample of ~ 10 4 early- type galaxies from the SDSS 
database. The SDSS sample includes photometric mea- 
surements in the u*, <?*, r*, i* and z* bands, as well as 
spectroscopic information. (See York et al. 2000 for a 
technical summary of the SDSS project; Stoughton et al. 
2002 for a description of the Early Data Release; Gunn 
et al. 1998 for details about the camera; Fukugita et al. 
1996, Hogg et al. 2001 and Smith et al. 2002 for details 
of the photometric system and calibration; Lupton et al. 
2001 for a discussion of the photometric data reduction 
pipeline; Pier et al. 2002 for the astrometric calibrations; 
Blanton et al. 2002 and Strauss et al. 2002 for details of 
the tiling algorithm and target selection.) 

As Bernardi et al. (2003a) discuss, the luminosities L, 
half light radii R and internal velocity dispersions a of 
the galaxies can all be reliably estimated from the data. 
Therefore, the sample is well-suited to estimating the dis- 
tribution of early-type galaxy stellar velocity dispersions. 
This is done in Section 2. Section 3 provides a simple es- 
timate of how the shape of the velocity function will be 
modified if the distribution of circular velocities from later 
type galaxies is added. Since we do not have measured 
circular velocities of late-type galaxies, the results of this 
section are, perhaps, less secure. Section 4 summarizes our 
results. 

2. THE DISTRIBUTION OF EARLY-TYPE GALAXY 
VELOCITY DISPERSIONS 

The SDSS main galaxy sample from which the Bernardi 
et al. (2003) early-type sample was compiled is magni- 
tude limited at both faint and bright apparent magni- 
tudes. In the r* band which we will use in this paper, 
14.5 < m r * < 17.77. Therefore, the observed distribu- 
tion of apparent magnitudes, sizes and velocities is not 
a fair estimate of the intrinsic distributions. Bernardi 
et al. (2003b, c) present a maximum likelihood analysis 
which accounts for the selection effects and measurement 
errors, at the cost of assuming a parametric form for the 
intrinsic distribution of luminosities L, sizes R a and veloc- 
ity dispersions a. In particular, Bernardi et al. assumed 
that the joint distribution of M = — 2.51og 10 L, log 10 R 
and log 10 a is tri-variate Gaussian. With this parame- 
terization, their analysis showed that the distribution of 
V = log 10 a is Gaussian with mean 2.20 and rms 0.11. 

Here we have chosen instead to show a nonparametric 
1/V m az estimate (Schmidt 1968) of the velocity dispersion 
function in Bernardi et al.'s sample. To calculate V max 
we assumed a fiat cosmological model with Q — 0.3 and 
Hubble constant of 70 km s _1 Mpc _1 , and used the SDSS 
apparent magnitude limits (14.5 < m. r » < 17.77). The 
lines with error bars in Figure 1 show this estimate. The 
Gaussian found by Bernardi et al. provides a reasonable 
but not perfect fit, so we have not shown it. As we discuss 
below, this is because they assumed that that the scatter 
around the mean a — L relation was the same for all L, 
whereas, in fact, the scatter depends weakly on L. The 
three curves in the top panel show different estimates of 
(j)(V) = ln(10) <j(f)(cr). They were obtained as follows. 
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Fig. 1. — Velocity dispersion function in the SDSS. Error bars 
show </>(c) for the early-type galaxy sample, estimated using the 
1 /Vmax method. Short dashed line which drops most sharply shows 
the result of transforming luminosities to velocity dispersions using 
the mean <r — L relation shown in the middle panel of Figure 2. 
Bold solid line shows the result of accounting for the scatter around 
the mean relation; it is well described by the form given by equa- 
tion (4), with parameters indicated in the panel. Long dashed line 
which provides a better fit to the data shows the result of convolving 
the intrinsic distribution shown by the solid line with measurement 
errors. Solid curve in bottom panel shows the ratio of the solid line 
in the top panel to the fitting formula, and error bars in bottom 
panel show log 10 of the ratio of the curve traced out by the error 
bars to the long dashed line in the top panel. 

The velocity dispersion function is 

cf>(a) = J dL<p(L)p(cr\L), (1) 

where p(a\L) is the distribution of a at fixed L. The joint 
distribution of a and L in this sample is shown in Fig- 
ure 2. In each panel, the contours show lines of equal 
probability, with levels chosen to be 1/2, 1/4, 1/8 etc. the 
value at the maximum. In the panel on the left, each 
galaxy was weighted by l/V max ; in this case, the joint dis- 
tribution p(a, L), is well described by a series of concentric 
ellipses whose major axes are all reasonably well aligned 
with each other. Notice that the probability distribution 
falls-off smoothly as one moves upwards and to the right 
along the major axis of the ellipse — there is no evidence 
of a sharp cut-off. The two dashed lines show Bernardi et 
al.'s maximum likelihood estimate of the two bisector fits: 
the average of the a — L and L — a regressions given in 
equation (2) below. The figure suggests that the distribu- 
tion of velocities at fixed luminosity is sligthly narrower 
for the brightest galaxies, a fact we return to later. 

Since we are more interested in p(p\L), the distribution 
of a at fixed L, than in the joint distribution p(a,L), we 
weighted each galaxy by l/[4>(L)V m ax], where 4>(L) is the 
value of <f> when that galaxy's luminosity is inserted into 
the luminosity function, and then re-plotted the joint dis- 
tribution of a and L. The middle panel of Figure 2 shows 
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Fig. 2. — Joint distribution of luminosity and velocity dispersion. 
In each panel, contours show lines of equal probability chosen to 
be 1/2, 1/4, 1/8 etc. the value at the maximum. In the panel on 
the left, each early- type galaxy was weighted by 1/Vmax, whereas 
the weighting was l/[tj> e (L) V m ax] and l/[(j> e (a) Vmax] in the middle 
and rightmost panels. The dashed line in each panel, and the text 
in the bottom, show Bernardi et al.'s (2003b) maximum likelihood 
estimates of these relations (sometimes called the bisector, inverse 
and direct fits). 



the resulting contour plot. The panel on the right shows 
the result of weighting each galaxy by 1 / [<p{cr)V m ax]', i-e., 
this panel shows p(L\cr). The dashed lines show 



log 10 er M„j = 2.2 - 0.102 (M„ + 21.15 + 0.85 z) 
Mr* log 10 cr\ =-21.15- 0.85 z- 5.75 (log 10 cr - 2.2$) 



from Bernardi et al. (2003b), where z is the redshift. (The 
values —21.15 and 2.2 are what Bernardi et al. estimate 
the mean values of M r * and log 10 a are in their sample.) 
The first of these relations is the inverse relation (a as a 
function of L, e.g. Schechter 1980) and is shown as the 
dashed line in the middle panel, whereas the second is 
the direct relation (L as a function of a), and is shown 
in the panel on the right. The lines in the panel on the 
left show the bisector fits, obtained by averaging the direct 
and inverse relations. (If the first of the relations above is 
V = (M — ainv)/binv and the second is M = ad ir + bdi r V, 
then the two average values of the slope are M <x V (bdir + 
b inv )/2 and V oc M {l/b dir + l/b inv )/2.) The text in the 
bottom of each panel gives the slopes of these relations. 
The fact that these slopes are rather different from one- 
another is a consequence of the fact that the a — L relation 
is not very tight. 

When luminosities are available but velocity dispersions 
are not, equation (1) shows that </>(cr) can be approximated 
if one has a good model of p[a\V). A first approximation, 
then is to assume that p(a\L) is sharply peaked about a 
characteristic value. Because we are interested in a at fixed 
L, this characteristic value is given by the first, rather than 
the second, of equations (2). Comparison of the middle 
panel of Figure 2 with those on either side of it, shows that 
this choice is important: the relations shown in the differ- 
ent panels are quite different from each other, so changing 
the relation used to transform from L to a will change the 



predicted p(a). Equation (1) shows that it is the relation 
in the middle panel which should be used. 

The short dashed curve in Figure 1 shows the result 
of transforming all observed luminosities into velocity dis- 
persions log 10 crpj using the first of equations (2) and then 
making a 1/V max estimate of the resulting velocity dis- 
persion function. Clearly, this procedure underestimates 
the number of objects with log 10 a > 2.3 by a large fac- 
tor. This shows explicitly that simply changing variables 
is an unreliable way to proceed. We show below that this 
happens because the scatter around the mean (a\L) is sub- 
stantial (see Figure 2). 

To incorporate the effects of this scatter, we added a 
Gaussian random variate with 

rms(F|M) = 0.079 1 + 0.17(M r „ + 21.15 + 0.85 z) (3) 

to log 10 <tfj for each galaxy. The solid line in Figure 1 
shows the resulting 1/V m ax estimate of <^(<r). It is in much 
better agreement with the actual distribution. The small 
discrepancy which remains is primarily due to observa- 
tional errors. To check this, we added a second Gaussian 
random variate to each log 10 crpj estimate, with rms given 
by the observational error quoted in Table 2 of Bernardi 
et al. (2003a). The associated distribution is shown by 
the long dashed line in the top panel of Figure 1. The er- 
ror bars in the bottom panel show the difference between 
the actual measured values (indicated by the error bars in 
the top panel) and the long dashed line in the top panel. 
This indicates that our method accounts quite accurately 
for the effects of intrinsic scatter and measurement errors. 
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Fig. 3. — Distribution of p(V\M) for a number of small bins 
in M . Error bars show this distribution when the actual velocity 
dispersions are used, and lines show the result of transforming the 
luminosities using the mean a — L relation and then accounting 
for scatter around the mean relation and measurement errors. A 
number of realizations of the scatter and errors are shown. Bottom 
panel shows the rms values; filled circles show the total observed 
scatter, open circles show the contribution from measurement errors, 
and crosses show the total scatter in each of the realizations. 
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The solid line in the bottom panel shows log 10 of the ratio 
of the solid in the top panel to the smooth fitting function 
described in equation (4) below. This indicates that our fit 
is an accurate description of the intrinsic cf>(a) distribution. 

The weak link in the procedure above is the assump- 
tion that the distribution around the mean a — L rela- 
tion is Gaussian (although the analysis in Bernardi et al. 
(2003b) suggests this should be reasonably accurate) with 
rms given by equation 3. To check this, the error bars in 
Figure 3 show ln(10) ap(a\L) for a few narrow bins in lu- 
minosity. The lines show the same quantity when log 10 ctfj 
plus intrinsic scatter plus measurement error are used in- 
stead of the actual log 10 a. We have actually shown several 
realizations of this procedure, so as to give some indication 
of how well 4>(a) can be determined from ~ 9000 galaxies. 
The filled circles in the bottom panel show the measured 
value of the rms scatter for each bin in L, and the open 
circles show the typical contribution from measurement 
error (the location along the x axis is given by the mean 
value of log 10 a in each bin). The crosses show the es- 
timated rms from each of the mock realizations. Except 
for the least luminous galaxies, our procedure recovers the 
observed p{a\L) distribution quite well. 

2.1. A fitting formula for 4>{°~) 

Although the velocity function can be described by 
transforming the luminosity function using the a — L rela- 
tion and its scatter, it is useful to describe 4>(a) in a way 
that is independent of waveband. The previous section 
showed that we had an accurate method for generating 
several realizations of the distribution of velocities. There- 
fore, we made several such realizations, measured <j){o~) in 
each, and then found the parameters (a, (3, er») which pro- 
vided the best fit to 



\a) da = <t>* ( — 



exp 



where 0» is the number density of galaxies. For this sam- 
ple, Bernardi et al. (2003b) estimate 0* = 0.002 (/i^Mpc)" 3 , 
whereas the 1/V max estimate, 0.0022 (/i^Mpc)" 3 , is slightly 
higher. We chose this functional form because, following 
Schechter (1976), the luminosity function is usually fit to 
such a function, but with (3 set equal to unity. The result 
of changing variables using (L/L*) = (ct/ct*)^ would re- 
quire 0(<t) = [^/T(a'/f3')}(a/a,) a ' exp[-(a/a*f} (3'/a 
which is of the form given above. Notice that, in this case, 
a' = atf) and f3' — flip: the parameters which specify 4>(a) 
are determined by fits to the luminosity function and to 
the a — L relation (e.g., Figure 2 shows that ip = 3.9). 
Since we have already shown that simply changing vari- 
ables is inaccurate, we chose to keep this functional form, 
but allowed a and f3 to be determined by the fit to 4>{a), 
rather than by the fits to (j>(L) and the a — L relation. 

The best fit values, (<f>, , a* , a, 0) = (0.0020±0.0001, 88.8± 
17.7, 6.5 ±1.0, 1.93 ±0.22), with </>» in (/i^Mpc)" 3 and a* 
in km s _1 , are indicated in the top panel of Figure 1. (Note 
that T(6. 5/1.93) w 2.88.) The value of logger* = 1.95 is, 
apparently, substantially smaller than the value 2.2 esti- 
mated by Bernardi et al. (2003b). This apparent dis- 
crepancy is resolved by noting that the mean value of a 
computed from equation (4) is ct* T[(a + 1)/ (3}/Y{a/ (3) « 
160 km s _1 . This is in excellent agreement with the mean 



value 10 2 ' 2 = 160 km s _1 derived by Bernardi et al. The 
solid line in the bottom panel indicates that this functional 
form describes the intrinsic 4>(a) distribution rather well. 

Although we have quoted rms values around the best-fits 
given above, the best-fit values are, in fact, strongly corre- 
lated with one another, suggesting that the data strongly 
constrain some combination of these parameters. Natural 
choices of such combinations are the mean, a, and most 
probable, a, values of a. Equation 4 shows that er* = 
aT(a//3)/T[(a + l)/0\ and (a-l)/(3 = (a/a*)?, which al- 
low us to provide a good approximation to the covariances 
between cr*, a and (3. This is illustrated in Figure 4: sym- 
bols show best-fit parameter values; symbol size indicates 
the goodness-of-fit. The smooth curves in each panel show 
(3 » (14.75/a) - 8 , and cr* = 161 T{a/(3)/T[{a + The 
fits return a = 160 ± 1.6 km s _1 values which only differ 
by about one percent, and the ratio a/a — 0.95 ± 0.007 is 
also very well determined. 

Although we do not show it, the scatter in is slightly 
correlated with a. This is because if one changes the shape 
of the probability distribution by allowing more galaxies in 
the faint end (i.e., by changing a), one must decrease the 
area under the curve at the bright end (because the area 
must integrate to unity) . Since it is the bright end which is 
better measured, it cannot change too much, which means 
</>* must increase to compensate. 

We have also performed the integral in equation (1) 
numerically, using the z = luminosity function from 
Bernardi et al. 2003b, and using equations 2 and 3 for 
p(a\L). Fitting equation (4) to the result gives a = 8 with 
cr*, and (3 as shown by the lines in Figure 4. 

The result of convolving the intrinsic distribution with 




Fig. 4. — Covariances between the best-fit parameters obtained 
in different Monte-Carlo realizations of the velocity distribution; 
symbol size indicates goodness-of-fit. The loeii traced out by the 
symbols in each panel are well described by requiring that all fits 
yield the same values for the mean a, and the most probable velocity 
dispersion a. Large filled circles show the best-fitting values we 
report in the previous figures. 
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measurement errors, modelled as a Gaussian distribu- 
tion with rms 0.035 in log 10 a, is well described by equa- 
tion (4) with (<7*,a,/3) = (88.8,6.5,1.8). The error bars 
in the bottom panel of Figure 1 show log 10 of the ra- 
tio of the observed 4>{a) to this form; clearly, this con- 
volved distribution provides a good description of the data. 
We conclude that our fitting formula, equation (4), with 
(cr*,a,/3) — (88.8,6.5,1.93), provides a good description 
of the intrinsic 4>{<r) distribution. 

3. THE CONTRIBUTION FROM LATE-TYPE GALAXIES 

So far we have studied the distribution of early-type 
galaxy velocity dispersions. The circular velocity v c is the 
analogous measure of the potential well of a spiral galaxy. 
Although we do not have measured values of v c for any 
of the SDSS galaxies which are not early-types, we can 
build a model of the contribution to the velocity function 
following the method used in the previous section: we first 
estimate the luminosity function of galaxies which are not 
early types, 4> ne (L), and we then use the v c — L relation to 
change variables from <j) ne (L) , being careful to account for 
inclination effects, which are expected to partially obscure 
the observed luminosities of late- type galaxies, and the 
intrinsic scatter around (v c \L). Finally, to compare with 
the velocity dispersion function of early-type galaxies, we 
convert from v c to velocity dispersion by assuming that 
a = v c / \/2. 

We estimated (j> ne (L) as follows. The luminosity func- 
tion of the entire SDSS sample, </> tot (L), has been esti- 
mated by Blanton et al. (2003). The SDSS photometric 
pipeline outputs a number of different estimates of the 
magnitude of a galaxy. Although Blanton et al. use Pet- 
rosian magnitudes, Bernardi et al.'s (2003b) estimate of 
the early type galaxy luminosity function did not. There- 
fore, we estimated the luminosity function of the Bernardi 
et al. early type sample using Petrosian magnitudes, 
<p e (L), and then set (j> ne (L) = (j> to t(L) - <p e {L)- (The Pet- 
rosian luminosity accounts for about 85% of the total lu- 
minosity in a deVaucouleurs profile. We found that <fi e (L) 
was indeed well approximated by simply rescaling all lu- 
minosities in Bernardi et al.'s 4>{L) by this factor.) 

Our next problem is to estimate the v c — L relation for 
later-type galaxies. Giovanelli et al. (1997) report that an 
inverse fit to the v c — L relation yields log 10 2v c = 2.5 — 
21.1/7.95- {Mi - 5 log 10 foioo)/7.94. Applying their fitting 
procedure to the early type galaxy sample we used in the 
previous section yields coefficients which are not formally 
the same as those of the (cr\L) relation of equation (2). 
However, the numerical value of the slope only changes 
from —0.102 to —0.104. Because these two are very close, 
we will assume that Giovanelli et al.'s fit approximates 
(v c \L) closely. 

Tully et al. (1998) also report a fit to the v c — L re- 
lation, and we have checked that it is very close to the 
one from Giovanelli et al. which we have chosen to use. 
We made this choice because Giovanelli et al.'s fit comes 
with a model for the scatter around the mean relation: at 
fixed velocity dispersion, the intrinsic scatter around their 
bivariate fit is ei nt = 0.26 — 0.28(log 10 2v c — 2.5) mags. We 
converted this into a scatter in log 10 a by dividing ej„t by 
7.94. Notice that this makes the scatter around the mean 
v c —L relation substantially smaller than it is around (c|L). 

To use these results, we first converted our simulated 



distribution of r* magnitudes into M° by setting Mf 
\ I ... — 0.9 (this is motivated by Fukugita et al. 1995 who 
suggest that the conversion factor is 0.95 for ellipticals, 
0.86 for SOs and 0.89 for Sabs). We then corrected lumi- 
nosities to face-on values following the discussion in Tully 
et al. (1998) (also see Giovanelli et al. 1995). This correc- 
tion makes use of the observed axis ratio b/a: namely, 

Mj = [Mf s + 3(16.9 + 5 log 10 MIA 1 - 5), where 
fl = -0.201og 10 (6/o). (5) 

In practice, galaxies are observed to have a range of axis 
ratios. Khairul-Alam & Ryden (2002) provide estimates 
of this distribution for SDSS galaxies, but we chose not to 
use their results because they do not account for selection 
effects. Our concern is prompted by the fact that estimates 
of p{b/a) in the Bernardi et al. sample, with and without 
l/Vmoz weighting, do differ from each other (see Figure 5). 
If the intrinsic axis ratio is tq , then 



p{b/a) = {b/a)J J 2 , (6) 
V 0/a) 2 - ?o 

if the distribution of inclination angles is random. This is 
the distribution of b/a values we chose to use. 

We made mock realizations of the contribution to the 
velocity function from objects which are not early-types 
by assuming that all galaxies which are not early-types 
have ro = 0.2 and follow the v c — L relation above. (In 
practice, there will be a range of tq values which our pro- 
cedure ignores, but the results to follow do not depend 
strongly on the precise value of ro.) The histogram which 
has the fewest galaxies with large values of a in Figure 6 
shows this estimate. Since the scatter around the mean 
v c — L is relatively small, accounting for it is not as im- 
portant as it was for the early- type galaxies. Similarly, 
correcting to face-on values by assuming that all galaxies 
have the same b/a — {b/a) (i.e., setting g = 0.056) and 
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Fig. 5. — Observed distribution of early-type galaxy axis ratios 
b/a (dashed line), and after weighting each galaxy by 1/V ma x- ac- 
counting for selection effects is important. 
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Fig. 6. — Estimate of the contribution to the distribution of 
velocities from late- type galaxies. The two histograms show the 
effect of transforming <p ne (L) and (frtot(L) using the Tully-Fisher re- 
lation and accounting for scatter around the mean relation as well 
as correcting for inclination and intrinsic absorption effects. The 
dashed line shows the result of transforming <j>tot(L) by ignoring 
the scatter around the Tully-Fisher relation and around the typical 
inclination and absorption correction. The solid line shows the con- 
tribution to the statistic from early- type galaxies; they dominate at 
a > 200km s" 1 . 

ignoring the scatter only results in a small underestimate 
of the distribution of large v c systems. This suggests that, 
for later-type galaxies, simply changing variables from L 
to v c should be reasonably accurate, provided one first cor- 
rects all luminosities to face-on values (but see discussion 
below). 

The solid curve which extends furthest to the right in 
Figure 6 shows the contribution from early type galaxies 
we discussed in the previous section; clearly, they domi- 
nate the statistic at a > 200 km s . To make this point 
even more clearly, the dashed line shows the result of as- 
suming that all galaxies are spirals, and so changing vari- 
ables in <ptot{L) rather than (j> ne (L), and the histogram 
shows the effect of including the scatter around the mean 
Tully-Fisher and inclination/absoprtion corrections. Our 
conclusion that early-types dominate at large velocity dis- 
persions is still valid. 

The correction to face-on values is large, on the order of 
0.5 mags, and, while well-defined, rather uncertain. This is 
why we have not performed fits to the contribution from 
later-types, nor have we tried to fit equation (4) to the 
sum of the two contributions (shown as the solid lines in 
Figure 6). 

4. DISCUSSION 

We have presented estimates of the distribution of ve- 
locity dispersions <P(<j) of early-type galaxies, and have 
shown that estimates which use the mean a at fixed L, 
(cr\L), to change variables from L to a, and ignore the 
scatter around the mean a — L relation underestimate the 



true number density of large velocity dispersion systems 
by large factors (Figure 1). We have shown that the de- 
pendence of virial velocity dispersion on galaxy luminosity 
is a power law (Figure 2), and we have derived an accurate 
model for the scatter in a at fixed luminosity (equation 3 
and Figure 3). Finally, we have provided a simple fitting 
formula for (f>(cr) (equation 4). 

We have also built a simple model of the contribution 
to the velocity function from galaxies which are not early- 
types. Our results suggest that, at velocity dispersions 
above about 200 km s _1 , early-type galaxies dominate 
the statistics (Figure 6). Thus, we have demonstrated 
that, at large a the velocity dispersion function falls as 
exp[-(cr/88.8 km s -1 ) 193 ]. 

The method we used for using observables other than a 
to estimate <j>(cr) is general. For instance, if one wishes to 
use the Fundamental Plane relation to derive <j>(<r) from 
photometric data only, then one requires knowledge of the 
mean velocity dispersion at fixed size and surface bright- 
ness, (cr\R , Ho), as well as the scatter around this mean 
relation. Note that it would be incorrect to use the coef- 
ficients of the usual direct fit to the Fundamental Plane 
relation, (i£ |er, (reported, e.g., in Table 2 of Bernardi 
et al. 2003c) to make the change of variables, for the same 
reason that it would have been incorrect to use the co- 
efficients of the (L\a) relation, rather than those of the 
(<t\L) relation (although, because the Fundamental Plane 
is tighter, the difference between the slopes will be smaller, 
and the effect of the scatter less pronounced.) 

Recent work (e.g. Trujillo et al. 2001) has revived in- 
terest in Sersic's (1968) generalization of the deVaucouleur 
(l/4)-profile to (l/n)-profiles. In particular, n appears to 
be rather tightly correlated with a. Use of this correla- 
tion may be a more promising way (than transforming the 
luminosity function) to estimate <f>(<r) from photometric 
information, and is the subject of work in progress. 

The mass function for clusters of galaxies cuts off sharply 
at large masses, as does the galaxy luminosity function. 
Our results indicate that we can now add the velocity dis- 
persion function to this list — a simple power law cannot 
describe the shape of (j>(a). As Schechtcr (2002) discusses, 
the dearth of galaxies with large velocity dispsersions con- 
tains important information about the gastrophysics of 
how the most massive galaxies must have formed. Indeed, 
Loeb & Peebles (2002) have used our measurement of </>(cr), 
in particular, our finding that values of a > 350 km s _1 
are extremely uncommon, to argue that the stars in the 
galaxies with largest velocity dispersions must have formed 
at sufficiently high redshift that gas dissipation effects are 
small. Kochanek (2001) has pointed out that combining a 
lensing based estimate of 0(c) with the one based on the 
motions of the stars, such as that presented here, provides 
powerful constraints on models of galaxy formation. By 
the time the SDSS survey is complete, a lensing based esti- 
mate of the velocity dispersion function should be possible. 
This will almost certainly measure velocity dispersions on 
larger scales than the few kpc scale probed by our mea- 
surement. Therefore, a comparison of the two will provide 
information about the effects of dissipation and baryonic 
contraction. 
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